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, Wilson loops in a turbulent fluid are shown to respect a speciflc area law cor- 

I responding to the Kolmogorov scaling. This law leads to the condensation of a 

• complex-valued scalar field minimally coupled to the velocity field. We use this find- 

ing to estimate a v.e.v. of the dual Higgs field, which appears in the hydrodynamic 
description of a spatially confining dual Landau-Ginzburg theory. The temperature 

■ dependence of all other parameters of this theory is found upon a comparison with 
the spatial string tension and the chromo-magnetic vacuum correlation length of the 

' Yang-Mills gluon plasma. In particular, a nonperturbative contribution to the shear 

' viscosity of the dual fiuid comes out exponentially suppressed with temperature. 

>: 

■ Interactions of the dual Abrikosov vortices with excitations of the fiuid yield a long- 

in ■ 

, range Aharonov-Bolim effect. This effect is shown to take place for all but calculated 

CN " discrete values of the product of the kinematic viscosity of the fluid to the coupling 

O! 

■ constant of the dual Higgs and velocity fields. In the presence of a Chern-Simons 
^ ■ term with the parameter 0, these discrete values remain the same for sufficiently 

^ I large values of 0. Furthermore, the Chern-Simons term leads to the appearance of 

■ - - ■ knotted dual Abrikosov vortices. 

I. INTRODUCTION 

The classic 't Hooft-Mandelstam scenario of confinement is based on the condensation of 
magnetic monopoles This phenomenon is described in terms of a magnetically-charged 
Higgs field. In case of the simplest, SU(2), Yang-Mills theory, the resulting Abelian Higgs 
model is just a 4-d generalization of the Landau-Ginzburg type-II superconductor, whose 
gauge field is dual to the A^-gluon of the initial Yang-Mills theory. It is also generally known 
that, in the deconfinement phase of the Yang-Mills theory, large spatial Wilson loops keep 



exhibiting an area-law behavior, which is the essence of spatial confinement (see Ref. Q for 
the corresponding lattice results in the SU(2) Yang-Mills theory). Analytically, spatial con- 
finement can with a good accuracy be described in terms of soft stochastic chromo-magnetic 
Yang-Mills fields jsl- Thus, above a certain temperature of the dimensional reduction, the 
soft part of the SU(2) Yang-Mills gluon plasma can be modeled by means of a dual Landau- 
Ginzburg theory with temperature-dependent couplings. 

On the other hand, it is a well-known property of a conducting turbulent medium that it 
expels an external magnetic field ^. Superimposing this fact with the 't Hooft-Mandelstam 
scenario of spatial confinement, one can guess that the dual-Higgs condensate, which within 
this scenario expels external chromo-electric fields, is formed owing to some dual turbulent 
medium. Such a medium implies an effective hydrodynamic description of the said dual 
Landau-Ginzburg theory, i.e. a description where the dual gauge field is substituted by 
a velocity field v^. Recently j^, the ratio of the shear viscosity to the entropy density in 
a medium of magnetically-charged hydrodynamic fiuctuations has been found to be much 
smaller than unity. We note also that a formally similar model, where the velocity plays a 
role of the gauge field interacting with quarks, was studied recently in Ref. j^. 

In this paper, we consider condensation of the dual Higgs field interacting with the dual 
velocity similarly to the formation of a heavy-quark condensate in QCD, which occurs due 
to the interaction of this quark with soft stochastic Yang-Mills fields [7|. To this end, we first 
find an area- dependence of the Wilson loop corresponding to the Kolmogorov scaling law in 
turbulence Q|. We consider then a complex- valued scalar field, which interacts via a covariant 
derivative with the velocity field. Using in the effective action of this complex-valued field 
the so-obtained Wilson loop, we get an estimate for the v.e.v. of the dual Higgs field. Using 
furthermore the known temperature behaviors of the spatial string tension and the chromo- 
magnetic vacuum correlation length in the Yang-Mills theory, we evaluate also the kinematic 
viscosity of the dual fiuid and the phenomenological coupling of the dual Higgs and velocity 
fields. This procedure fixes completely the temperature dependence of all the parameters 
entering the effective hydrodynamic description the dual Landau-Ginzburg theory. Having 
these parameters fixed, we proceed to the analysis of topological effects stemming from 
the interaction of excitations of the dual fluid with the dual (i.e. carrying electric fluxes) 
Abrikosov vortices j^, which are present in the dual Landau-Ginzburg theory. To this end, 
we perform a duality transformation of the Wilson loop describing such an excitation, and 
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find in particular a long-range Aharonov-Bohm-type interaction between excitations of the 
dual fluid and the dual Abrikosov vortices. Finally, we extend this analysis to the case where 
a Chern-Simons (CS) term is present. 

The paper is organized as follows. In the next Section, we derive for the Wilson loop in 
hydrodynamics a novel area law corresponding to the Kolmogorov scaling. In Section III, 
we consider an interaction of the velocity field, which respects this scaling, with a complex- 
valued scalar field. Using the obtained specific area law for the Wilson loop, we evaluate 
a v.e.v. of the dual Higgs field in terms of the energy dissipation rate and the kinematic 
viscosity of the dual fiuid. We make further a natural assumption that the dual Landau- 
Ginzburg theory represents a dimensionally-reduced high-tcmpcraturc phase of a 4-d dual 
Abelian Higgs model. In the London limit, this model yields a string tension and a vacuum 
correlation length, which can be compared with the known spatial string tension and the 
chromo-magnetic vacuum correlation length of the Yang-Mills gluon plasma. Using this 
correspondence, we obtain the kinematic viscosity of the dual fiuid in terms of the following 
three quantities: the energy dissipation rate, the temperature, and the finite-temperature 
Yang-Mills coupling. In the same way, we also obtain in terms of these three quantities a 
dimensionful phenomenological coupling /x of the dual Higgs and velocity fields. In particu- 
lar, the obtained purely nonperturbative kinematic viscosity turns out to be exponentially 
suppressed with temperature. 

As was mentioned above, the Landau-Ginzburg theory at issue possesses dual Abrikosov 
vortices, which interact with excitations of the dual fiuid. In Section IV, we show that 
among these interactions there exists a long-range Aharonov-Bohm one, and find a relation 
between the kinematic viscosity and the coupling fi, for which this interaction disappears. 
In Section V, we additionally consider the efi^ects produced by the CS term. First, we 
briefiy show that, in the absence of the dual Higgs field, the CS term yields a self-linking 
of the contour along which the dual excitation evolves. Then we consider an excitation 
of the dual fiuid in the full theory, which includes also the dual Higgs field, and obtain 
there the Aharonov-Bohm-type interactions in the presence of a CS term. In particular, 
at sufficiently large values of the ©-parameter entering the CS term, we obtain an analytic 
expression for the Wilson loop associated with an excitation of the dual fluid. Using this 
expression, we find a relation between the kinematic viscosity and the coupling /i, for which 
the Aharonov-Bohm effect in the large-© limit disappears. Remarkably, this relation turns 
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out to be the same as the one obtained in Section IV in the absence of the CS term. We 
also notice that the temperature behaviors of the kinematic viscosity and the couphng fi, 
obtained in Section III, allow this relation to hold, that makes it nontrivial. Furthermore, 
we explicitly find the appearance of knotted dual Abrikosov vortices in the same large-0 
limit. In Section VI, the summary of the results obtained is presented. In Appendices A 
and B, we provide some technical details of the calculations performed. 



II. A NOVEL AREA LAW OF THE WILSON LOOP IN TURBULENCE FROM 

THE KOLMOGOROV SCALING 

We start with considering a Wilson loop in turbulence {lo| . 

{W{C)) = (exp (^^ J d'xv,j,^ ^ , (1) 

and derive its area dependence corresponding to the Kolmogorov scaling law. In Eq. ([1]), 
j^(x; C) = §^ (ix^(r) 5(x— x(r)) is a conserved current, and v is the kinematic viscosity of the 
fluid, equal to the ratio of the shear viscosity to the density. The Kolmogorov hypothesis 8| 
states that the mean dissipation of the kinetic energy by a unit mass of the fluid, occurring 
in the unit time, is constant: 8 = const. Since S = v'^/t, where t ~ one gets the 

Kolmogorov scaling law: 

V ~ (^|x|)v^ (2) 

where v = |v|. This law can be given a field-theoretical description by evaluating v as 
^(f^)jr. Here {■ ■ effectively substitutes the average (■ ■ ■ ) from Eq. ([1]) at sufficiently 
small distances, where the Kolmogorov scaling holds. At such distances, one can neglect 
connected velocities' correlation functions higher than the two-point one, and approximate 

(WiC))^ as 

{W{C))^ ~ exp £dx, £dy,{v,{^)vM)T^ = exp (^-const " ^ 5^^') , (3) 

where S is the area of the minimal surface encircled by the non-selfintersecting contour C. 

The average (■ ■ ■ )jr can be formally defined by representing in the form = j^j^J-'uti, 
where J^,^n is a space-independent matrix with the Gaussian distribution: 

(■■■). = J^2[Il I ^^.^e-^l (■■■)• (4) 
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Here, in the general case of a (i-dimensional Euclidean space, n = {<P-d)/2, and the average 
is normalized by the standard condition = 1. Noticing that 



+00 -p2 



d 



n/2 



we can obtain the coefficient c in the formula {J^^yJ^xp) jr = c ■ {5^\5yp — 6^p6,^x). Indeed, this 
formula yields {J-''^^)jr = c ■ {d"^ — d), that leads to the value of c = ^^r^ = ^. Therefore, 
setting d = 3, one has 



(^J>^ = SI (-^-m-^Am)^ = 7|x 



2/3 



|x|4/3 

where the above-obtained mean value {J^^uJ^Xp) jr = |(5/^A<^!^p — ^np^ux) has been used at the 
final step. Equating now ^ {'^p}jr = 7^^^|x|-'^/'^ to the right-hand side of Eq. ([2]), we obtain 
from the Kolmogorov scaling law the coefficient 7 in terms of the energy dissipation rate: 



7 ~ £2/3_ 

Furthermore, using the correlation function 

7 XxVp 



{Vf,{^)vM)T = 2 |x|2/3|y|2/3 ('^ApV - ^XuSpp) , (5) 



one gets for the Wilson loop ([T]) the following expression: 



{W{C))^ = exp ( ji^dx^ j^dy^ {vp{^)vy{y)) ^ 

xy jT dx^Xy jT dy^y^ 

2/3 



For a fiat contour, the expression in the square brackets can be estimated as [■ ■ ■] ~ S^^^. 
(In particular, an elementary calculation yields for a circle [■■■]= 47r^/^5'^/^.) Thus, expres- 
sion ([3]) for the Wilson loop is recovered. 



III. A HYDRODYNAMIC DESCRIPTION OF A SPATIALLY CONFINING 
DUAL LANDAU GINZBURG THEORY 

A. Condensation of a complex scalar field interacting with the turbulent fluid 

Let us consider a complex-valued scalar field of some mass M, interacting with the 
turbulent fluid. The action of such a field has the form S = J d^x + M^|0p), where 
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the covariant derivative reads = dfj_ — ^v^. In the course of its spatial evolution, the field 
(j) yields the Wilson loop ([1]). We show how the specific area law in turbulence, Eq. Q, leads 
to the condensation of the field 0, i.e. to the formation of the condensate (From now 

on, we denote by (■ ■ • ) the average (■ ■ ■ 

To this end, we integrate in the partition function over the (f)- and the (/)*-fields, that 
yields the one-loop effective action 

-Dl + M^l 



We use here the normalization condition (r[0]) = 0, which corresponds to our assumption 
that the 0-field is condensed because of its interaction with the velocities f^'s respecting the 
Kolmogorov scaling. One can further represent (r[f^]) as the following path integral: 



(4vrs)3/2 



(6) 



Here (yV{C)) is given by Eq. ([3]), the factor of 2 stems from the complex- valuedness of 0, 
and y is a 3-d volume occupied by the system. The condensate of the 0-field can then be 
found from the relation 

(I0r) = -^^(rk]). (7) 

We calculate the path integral that enters Eq. ([6]) by reducing it to the one in an auxiliary 
space-independent electromagnetic field of a geometric origin. To this end, we use the 
fact that, for a sufficiently small non-selfintersecting contour C at issue, the area S of the 



minimal surface encircled by such a contour can be represented in the form S = y S2^/2, 
where E^j, = §^dx^Xv is the so-called tensor area |lo|. Furthermore, we use an elementary 
saddle-point integral 

_ 2/3 /S" _a2_ 2k. 

e'' ~W— / dXe 33727 where /t > 0, 







V TT 

owing to which we can represent the Wilson loop (E]) in the form 



In this formula, we have denoted 



{W{C))c^\l- / ciAe-^'-^(^)-^-. (8) 







^ , , , / const \^^'^ 8 , . 
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where "const" is the same as in Eq. ([3]). Then the S^,^-part of the exponent in Eq. ([8]) can 
be represented as an integral over the said auxihary space-independent electromagnetic field 




and b"^ = I 



Furthermore, the mass parameter M 



-^(A)-z.- 2 III/ dB I e~16G{A)-2 

Accordingly, the path integral in Eq. ([6]) becomes reduced to the so-called Euler-Heisenberg- 
Schwinger Lagrangian for a scalar particle jUl] (for reviews, see 12|): 

(47rs)3/2 sin(as) sinh(6s) ' ^ ^ 

x^{0)=x^{s)=0 

Here a and b can be expressed through the auxiliary electric and magnetic fields correspond- 
ing to the Abelian strength tensor as = | - + ^/(E^ - H^)^ + 4(E ■ H)^ 

H2 - E2 + ^(E2-H2)2 + 4(E-H)2 
should be sufficiently large in order for the mean size of the trajectory to be as small as 
(9(z/) to ensure the validity of the Gaussian approximation used in the derivation of Eq. Q. 
Therefore, we set M ~ proving afterwards that this parameter is indeed getting expo- 
nentially large with temperature. For such large masses in Eq. ([6]), one can safely expand 
the Euler-Heisenberg-Schwinger Lagrangian up to the leading nontrivial term as 

^1 + ^^52^. (11) 

sin(as) sinh(6s) 6 '"^ 

Then Eqs. ([6]), ([7]), ffTOj) . and (fTTj) yield for the condensate of the 0-field the following 
intermediate expression: 



The Gaussian 5^^-integration here can be readily performed as 

"+00 \ b2 foo 




dB^^ e'^m^ ■^Bl^ = A-K j dBB^e'^ = 192y27r=^/^[G(A)] 



5/2 



leading to a simple formula 

/o poo 

Recalling now the definition of G{X) given in Eq. (I9l), we obtain 

, _ const^/^ £ dX 
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In order to estimate the cut-off parameter Amin, we should go baclc to Eq. ([8]), where the 
A-integration first appeared. The saddle-point value of A, following from the integral in 
Eq. (ED, reads As.p. ~ ^ . Since the minimal size of the trajectory is estimated as z/, we 

conclude that As. p. > {8vY^^. Therefore, since the A-integral is "saturated" by the vicinity of 
As.p., we can set Amin = {SuY^^. Recalling that M ~ ^, we finally obtain for the condensate 
of the 0-field the following estimate: 

(|0p)oc^-ln-^. (12) 



B. A spatially confining dual Landau Ginzburg theory 

Let us now consider a dual Landau-Ginzburg theory, which can describe the phenomenon 
of spatial confinement. This theory contains a dual Higgs field and a dual gauge field. The 
dual Higgs field describes a 3-d medium of magnetic monopoles above the temperature 
of dimensional reduction. Such monopoles are therefore assumed to be those degrees of 
freedom whose condensation provides spatial confinement of gluons in the dimensionally- 
reduced high-temperature Yang-Mills plasma. 

We consider the simplest, SU(2), case, and adopt an effective hydrodynamic description 
of the corresponding dual Landau-Ginzburg theory. Within this description, the dual gauge 
field is replaced by the velocity field, which is further assumed to be represented as a sum 
of the low- and the high-energy modes. In a sense, such a separation of the velocity into 
two components is similar to that of the Landau model of superfiuid quantum liquids 13] . 



In our case, the high-energy modes have the standard, for the Landau-Ginzburg theory, 
Maxwell-type kinetic term F'^^[v\^ where -^^i^f^^] = d^v^ — d^v^. (Apparently, is only 

non-vanishing for a non-potential motion, for which the high-energy part of the velocity 
cannot be represented as a gradient of some scalar function.) To minimize the number 
of free parameters, we assume for this kinetic term a normalization factor of Rather, 
the average over the low-energy modes is given by Eqs. (jl]) and ([5]). It yields Kolmogorov 
turbulence for these modes, and can lead to the formation of the condensate ( IT2l) . 

We make now our main conjecture that the Kolmogorov turbulence of the low-energy ve- 
locity modes can be the origin of the dual-Higgs condensation in the dual Landau-Ginzburg 
theory. As a physical motivation for this conjecture serves the known fact that a conducting 
turbulent medium expels external magnetic fields j^. In the presence of a Higgs field, one 



9 



can introduce a Meissner-type description of this phenomenon by assuming that turbulence 
leads to the condensation of this Higgs field. Here we impose a dual version of this sce- 
nario, and use Eq. ( IT2l) as an estimate for the v.e.v. of the dual Higgs field in the dual 
Landau-Ginzburg theory: 

'73d ~ ^Js -In-^. (13) 

We introduce further a phenomenological coupling fi of the dual Higgs field and the high- 
energy velocity field v^. That is, the covariant derivative acting on the dual Higgs field 
reads Dx = dx — ifivx, and the coupling /i has the dimensionality of (mass). Furthermore, 
we consider the London limit of the dual superconductor. In this limit, the radial part of 
the dual Higgs field goes over to rj^d, and only the phase 6 of this field remains dynamical. 
Accordingly, the action of the dual Landau-Ginzburg theory in the London limit reads 

^3d = I d'x {^i^iM + + f^v.f^ . (14) 

In order to fix the temperature dependence of /i and u, we consider now a 4-d dual Abelian 
Higgs model, which yields the dual Landau-Ginzburg theory upon the high-temperature 
dimensional reduction. In the London limit of interest, the action of such a dual Abelian 
Higgs model reads 

Here is a dual gauge field, F^y[B\ = d^By — d^B^ is its strength tensor, is a dimen- 
sionless magnetic coupling, and we make no difference in notations between the 4-d and the 
3-d indices. The action of the corresponding dimensionally-reduced theory follows upon the 
replacement J d^x P j d^x, where /3 = 1/T. Matching now the fields and the couplings 
of the actions 84^ and S^^, we obtain the following three relations: 

Pvld = vld^ 9mBf, = fiv^, B^ = y^^M- (15) 

Let us now consider the soft part of the Yang-Mills gluon plasma which, at sufficiently 
high temperatures, also undergoes the dimensional reduction The two quantities 

characterizing this dimensionally-reduced soft part of the gluon plasma are the spatial string 
tension cr^ and the inverse spatial vacuum correlation length These quantities can 

be used to obtain the parametric temperature dependence of fi and u. To this end, we 
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notice that, since both (jg and characterize soft stochastic background fields, they are 
proportional to the corresponding power of g'^T, where g is the finite-temperature Yang-Mills 
coupling. Namely, one has jsl 

a, oc(/T)2, r'oc/T. (16) 

On the other hand, in the London limit of a dual superconductor, the string tension is given 
by the formula cr^ oc rjl^ In k, where k is the ratio of the dual-Higgs and the dual-vector-boson 
masses fl^ . We notice that, although the formal derivation of this expression implies the 
inequality In/t ^ 1, lattice simulations in the Abelian-projected Yang-Mills theory suggest 
the values of Iuk < 5 (cf. jl^ and references therein). Therefore, we have 

(Ts ~ ??4d = Tril^ ~ ■ In (17) 

where Eq. (fT3|) has been used at the final step. Using now Eq. (fT6|) . we obtain 

~ £ . In i- ~ g'T. (18) 

(We notice that, for presumably an accidental reason, the right-hand side of this relation 
coincides with the known ultrasoft scale, the inverse of which is the mean time needed for 
a parton undergoing Coulomb scatterings in the gluon plasma to defiect by an angle of the 
order of 1 15|.) Recalling now that £ = const in the initial Kolmogorov hypothesis, we 



obtain the sought temperature dependence of the dual kinematic viscosity: 

^ ^ _g-const.V. (19) 

Thus, we conclude that the soft contribution to the kinematic viscosity of the dual fiuid in our 
approach is suppressed with temperature as e"™"^*^* ". Consequently, the soft contribution 
to the shear viscosity of the dual fiuid is suppressed by the same factor. This finding 
parallels the smallness of nonperturbative contributions to the dual shear viscosity, which 
was found within other approaches [5]. We notice that, on the general grounds, perturbative 
contributions to the shear viscosity are also expected. Those should keep the ratio of the full 
shear viscosity to the entropy density above the lower l/(47rVbound which, in the isotropic 



case at issue, is predicted by the AdS/CFT correspondence [l6|. 

Comparing now the inverse vacuum correlation length in the dual Abelian Higgs model 
with that in the high-temperature Yang-Mills theory, we can get the temperature depen- 
dence of the phenomenological coupling /x. Indeed, the inverse vacuum correlation length in 
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the London hmit of the dual superconductor is given by the mass of the dual vector boson, i.e. 
~ OmV^d- The magnetic coupling can be obtained from the second and the third rela- 



tions ( fTSl) : gm = jJ^y/^- Using for ?74d expression ( fTTl) . we have ~ /iyf^z/ ■ In Further- 
more, since In^ ~ [cf. Eq. flTS]) ]. we can write the previous estimate as ~ g^fiy/Tu. 
Equating now this expression to the Yang-Mills one, Eq. f fT6l) . we get a simple formula 



/i~y^- (20) 

Using finally Eq. (fT9|) . we obtain the temperature dependence of /i in the form 

v^e™'^'^*-^, (21) 



where "const" is the same as in Eq. (fT9|) . 

Anticipating the results of the next Section, we will encounter there the relation ^ 
^ — , which should not hold in order for the Aharonov-Bohm-type interactions between 

integer ' 

the excitations in the dual fluid and the dual Abrikosov vortices to be nontrivial. From 
Eqs. ([IHD and we see that 

—const- 



/iz^ ~ Y e 2£ ^ 

i.e. it is an exponentially small quantity, which therefore may indeed be oc 1/integer. This 
makes the above-quoted constraint on the product fiu nontrivial. 



IV. THE AHARONOV BOHM EFFECT IN THE ABSENCE OF A 

CHERN SIMONS TERM 

A phase factor (i.e. an unaveraged Wilson loop) associated with an excitation of the dual 
fluid is defined in the same way as in Eq. ([1]), i.e. as exp(^ /x'^m^m)' where from now on we 
use short-hand notations J^ = J d?x and Jp = J (^^- Let us first consider the case without 
condensation of the dual Higgs field, i.e. set in Eq. (fl^ 773^ = 0. Then the phase factor 
averaged with the remaining Maxwell action, ^ f^F'^ul'^]' yields the following Wilson loop: 

<"''<^'> = ^''p (- 8^; i "'^ i ''^'R^) • (2^) 

We consider now the full action (fl^ . with 773^ 7^ 0, and average the phase factor 
exp(^ f^Vfj_j^^ with this action. That is, we calculate the following Wilson loop: 

{w{C))= [ vv^vevee--^^l^^''^^''^^'''^^^'^-^''''^^'->^^^l (23) 
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Here, 1] = risd, and the full phase 6 of the dual Higgs field is represented as a sum 6 = 6 + 6, 
with 6 experiencing jumps by 27i when going around dual Abrikosov vortices, while 6 being 
a Gaussian fluctuation around 6. The said jumps of 6 lead to the noncommutativity of two 
derivatives acting on this field: 

{d^d, - d,d^)6 = 27Te^,xJx. (24) 

Here J\ is the current of dual Abrikosov vortices, which populate the dual-Higgs vacuum. 
One can expect that the trajectory C, along which an excitation of the dual fluid evolves, 
links to the dual Abrikosov vortices. Such a linking should correspond to a long-range 
Aharonov-Bohm-type interaction between the said excitation and the dual Abrikosov vor- 
tices. 

To visualize such an Aharonov-Bohm interaction, we perform now a duality transforma- 
tion of the Wilson loop (123|) . To this end, it is first convenient to introduce two auxiliary 
fields as follows: 

The subsequent integration over 6 leads to the constraint d^Cf^ = 0, which can be resolved 
by representing as = s^^xdyifx, where the field y?^ has the dimensionality of (mass). 
Accordingly, = where = df^ip^ - d^ip^, and i C^,dJ) = 2-ni (ff,J^, where at 

the last step we have used Eq. Thus, the Wilson loop (123|) takes the form 

Notice that, throughout this paper, we work at the entirely classical level. For this reason, 
the Jacobian corresponding to the change of integration variables ^ — )■ is omitted, and the 
measure DJ^ in the functional integral has only a statistical (rather than a field-theoretical) 
meaning of counting vortices in their given configuration. 

Owing to the representation ( l25l) for the Wilson loop, the t;^-field becomes a 
Lagrange multiplier, which leads, upon the integration, to a functional 5-function 
^ {^^lu\^u{G\ + ^ipx) + ^i^l^. The G^-integration amounts to substituting G^, stemming 
from this (5-function, into -^x^l. Such a reads = —fif^j, — ^^^_lux fy ^^Dq^J^, where 
Dq^ = l/(47r|x — y|), j| = jx{y] G), and the conservation of has been used. Accordingly, 
the Wilson loop takes the form 

{W{G)) = I pj^p^^e^4-^*-+2-^^"^''-t('^^M+i^.../,9.^^o"..")T^ 
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or, equivalently, 
where 

K; ^ 27, j; + z^e^,, I dlDl^jl (26) 

The remaining (^^-integration amounts to introducing a rescaled field = V^m/('7"\/2), and 
denoting 

= 2^x^7fv. (27) 

That yields 

where Z)|^^ = e~'^'^^^'/(47r|x — y|). Thus, the Wilson loop becomes 

{W{C)) = e-^^->«^'2^'^o' j PJ^e"'''^->^-^^^'^^'. 

The last exponent in this formula can be simplified (cf. Appendix A), that yields the 
following result: 

{W{C)) = e-^^-,.^'^^''^-' J VJ^e-'^^^'^^'^^^y^^^-^'^'^-'+l^^''^^'^^^^^ (28) 

where L{j,J) = e^i^x J^J^O^Dq^ is the Gauss' linking number of the contour C and a 
dual Abrikosov vortex. The corresponding exponent e in Eq. (1281) describes a long- 
range Aharonov-Bohm-type interaction of the magnetic excitation with the dual Abrikosov 



2-iTi 



vortex. Thus, this interaction is nontrivial (i.e., e '^^ '^^"''"'^ ^ 1) as long as 

/xz/^— (29) 
mteger 

We notice that, at the formal level, Eq. (I28l) represents a 3-d counterpart of the corresponding 



4-d result for the Wilson loop in the Abelian Higgs model (cf. Ref. 17l|). 



V. THE AHARONOV BOHM EFFECT IN THE PRESENCE OF A 

CHERN SIMONS TERM 



We extend now the analysis of the previous Section to the case where the CS term is 
included. As a warm-up, we again start with the theory in which the condensation of the 
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dual Higgs field is absent, i.e. 77 = 0. The Wilson loop in such a theory is given by the 
formula 

{W{C)) = J Vv^e^ ■^^^i^^^"^"^^'^^'"'^''''^'"''^'^''''^'''^, 

where the dimensionality of the new parameter 6 is (mass)^. Imposing the gauge-fixing 
condition d^v^ = 0, we obtain the saddle-point equation 

—d'^v^ + ime^y\dyV\ = ij^, where m = 29z/. 

Seeking a solution to this saddle-point equation in the form = 1/^ + zV^, we get a system 
of equations 

d^U^ + me^^xduVx = 0, -d^V^ + me^.^d^U^ = j^. (30) 
The first of these equations can be solved with respect to f/^ as 

Ul = me,,, I Dl^dlVl (31) 



y 



Differentiating the second equation fl30p . and applying the maximum principle, one gets 
<9/iV^ = 0. Using this relation, one further obtains from Eq. (13T]) : e^uxd^Ux = mV^. Substi- 
tution of this formula into the second equation (130|) yields for that equation a remarkably 
simple form (—9^ + m^)V^ = j^. Therefore, one has = J^D^j^, while f/^, given by 
Eq. ([31]), can be calculated by virtue of Eq. (jM]), and reads = ^e^^^x jyiD^^ - D^)dyjl. 
Altogether, the resulting Wilson loop has the form 



(32) 

I x,y L '"■ 1 ) 

This expression generalizes Eq. (122]) to the case of O 7^ 0. Clearly, the G-term leads to the 
appearance of a self-linking of the contour C, as well as of a short-range self-interaction of 
this contour by means of the Yukawa propagator D^. We also notice that, when — in 
Eq. fl32]) . Eq. (122]) is recovered. Indeed, in this limit, one has ^(Dq^ — D^) — so that 

il/.3id:{D^o' - Dni) = il/^^o' - D^^odiji ^ l^l^r.difx = 0, 

since f jf, = 0. 

We proceed now to the duality transformation of the Wilson loop in a full theory, where 
the condensation of the dual Higgs field does take place, i.e. 1] 0. The corresponding 
generalization of Eq. ( 123]) reads 

{W{C)) = j Vv,V9V9 e~ •^''^^^^'•'^''^^''^^'^^'''^^^'''''^^^'^^'''''■^ (33) 
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The transformation leading from Eq. (123|) to Eq. (125|) remains the same, so that the coun- 
terpart of Eq. (l25l) in the presence of the CS term reads 

Because of the CS term, now ceases to be a Lagrange multipher. Nevertheless, since the 
^^-integration is Gaussian, it can be performed exactly. 

The corresponding saddle-point equation for reads Sfj^uxdyVx = ^k^, where we have 
denoted = e^y\dy{G\ + jJi^x) + ^j^. Owing to the conservation of fc^, a solution to this 
saddle-point equation reads v'^ = ^e^uxd^ ly^o^^x- Plugging this solution back into the 
exponent e*'^^; '^mC^m-Q^mi'A^i'I'a)^ a^j^fj using the above explicit expression for fc^, we obtain, upon 
some algebra, the following formula for the Wilson loop: 

X f VJ Vlo VG e-^==l~^'^^~8^*^''''"4^^''''^[*^''^''^'^^^^^'^^''^''^'^'''^"''^^]'^2i^ 

(34) 

Here, the argument of the first exponent coincides with the term containing the Gauss' self- 
linking number of the contour G, which was present already in Eq. (l32l) . In general, such 
self-linkings are always induced by the CS term 18|]. In addition, the functional integral in 



Eq. (iMj) describes interactions of a dual-fluid excitation, which evolves along the contour C, 
with the dual Abrikosov vortices, as well as self-interactions of this excitation in the presence 
of the CS term. 

In order to visualize all these interactions, let us perform the G^-integration first. Repre- 
senting the saddle-point expression for G^ in the form G^ = + i^^, we obtain a system 
of two saddle-point equations: 

e^yxdyLx - mN^ + = 0, e^^xduNx + niL^ = 0, 

where we have denoted = ^e^^xduVx + ^J/^- Owing to the conservation of n^, we find a 
solution to these equations in the form 

LI = -e,ux [ D^Jdlnl = m f D^^nl. 

Plugging the corresponding saddle-point expression for back into Eq. we obtain, 

after some algebra, the following general result: 
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For exponentially small z/'s of interest [cf. Eq. fll9p ]. the remaining (/^^-integration in 
Eq. (I34p should reproduce Eq. (1221) . which follows directly from the initial Eq. (I33p : 

{W{C)) j Vv^e--^-(^^^-'->''^^) = e-^^-.«^'^^'^«'. (36) 

Therefore, the (z/ — )■ 0)-limit can serve as a check for the above-performed calculation. The 
right-hand side of Eq. ( !35l) simplifies in this limit to the form 

and the Wilson loop becomes 

The Gaussian (/^^-integration in this formula yields, upon some algebra. 



Recalling the normalization of the integration measure "DJ^, discussed in Appendix A, we 
indeed recover the expected result (l36l) . Thus, our check of Eq. ( l35l) was successful. 
We consider now large values of the B-parameter, namely such that 

e > ^ = g^TS ■ e™'^^*-^^, (37) 

where Eq. (fT9l) has been used at the final step. For such O's, the action in the exponents 
on the right-hand side of Eq. (135|) becomes local, and the Wilson loop (134|) gets simplified 
to the form 

{W{C)) -> e~^~''-^''+3^^'^'^^(^-.«^"^''^^^o''"^^--''''^-^'^)x 



Furthermore, the (/^^-integration in this formula can also be performed analytically in the 
same limiting case ( 1371) . Referring the reader for the details to Appendix B, we present here 
the final result of this integration: 

{W{C)) e"^-^-^''"^^'^''^-^-^'''^''^'^x 
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X 



,2„2 



In this formula, Ai = ^^-q-, and i?^ = 27rJ^ + ^j^. Remarkably, in the limit (1571) . the 
initial CS term for the velocity, iQe^uxV^duVx, leads to the appearance of its counterpart 
■^^nuxiiidvix for the current j^, while the self-linking of the contour C, described by the 
first exponent in Eq. disappears. Rather, we observe the appearance of self-linkings of 
the dual Abrikosov vortices, as well as their linkings with the contour C, as described by the 
term ^^^e^y\J^ (^^a + m^V) ^u^o^ Lagrangian. In particular, its part ^^j^e^^xJ^ ■ 

£-jyg^]j^y yields the term — ^L{j, J) in the action. Thus, we conclude that, in the presence 



of the CS term, the Aharonov-Bohm-type interaction of the dual excitation with the dual 
Abrikosov vortex is nontrivial in the limit fl371) as long as 



/xz/^— (40) 
mteger 



Remarkably, this condition coincides with Eq. (l29l) . which was obtained in the absence of 
the CS term. In a similar way, the above term ^^^2 £^uX J\dl Dq^ means that the CS term 
makes dual Abrikosov vortices knotted unless ^ 7^ . f ^ . 

B ' integer 



VI. SUMMARY 



The 't Hooft-Mandelstam scenario models spatial confinement in the SU(2) gluon plasma 
by means of a dual Landau-Ginzburg theory. We have used an effective hydrodynamic 
description, where the dual-Higgs condensate is formed owing to a turbulent fluid, which 
respects the Kolmogorov scaling law. Such an approach yields an estimate for the v.e.v. of 
the dual Higgs fleld, Eq. ( fT3l) . Using this flnding along with the known expressions for the 
spatial string tension and the correlation length in the high-temperature Yang-Mills theory, 
we have found the parametric temperature-dependence of the dual kinematic viscosity v and 
the phenomenological coupling /i of the dual Higgs and velocity flelds, Eqs. (fT9|) and (12T|) . 
In particular, the product of these two quantities comes out exponentially suppressed with 
temperature, and hence parametrically smaller than unity. 

We have further considered a Wilson loop describing an excitation in the dual fluid. Such 
an excitation is involved, in particular, in a long-range Aharonov-Bohm-type interaction 
with the dual Abrikosov vortices, which are present in the dual Landau-Ginzburg theory. 
We have shown that this interaction is nontrivial for all but discrete values (|29|) of the 
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product /iz/. These values are smaller than unity, much as the calculated product /iz/ in 
general (cf. the previous paragraph), that makes constraint (129|) nontrivial. 

The Wilson loop describing an excitation of the dual fluid has further been calculated 
in the presence of a CS term. In the absence of the dual Higgs fleld, the result is given by 
Eq. ( l32l) . which includes self-linking of the contour of the Wilson loop. In the presence of the 
dual Higgs fleld, the result is given by Eqs. (l34l) and (l35l) . where the remaining y?^- integration 
can be performed analytically only in some limiting cases. In the case of a vanishingly small 
viscosity v at flxed B, the expected result (l36l) is reproduced. In another limiting case, 
when B's are as large as those given by Eq. ( 137|) . the y9^-integration can also be performed 
analytically, yielding Eq. (l39l) . This result contains terms, which describe self-linking of 
dual Abrikosov vortices, as well as their linking with the contour of the Wilson loop. The 
latter linking represents the same Aharonov-Bohm effect as in the absence of a CS term. In 
particular, the critical discrete values of the product /ii/, for which this effect disappears, are 
given by Eq. ( l40l) . that coincides with its counterpart ( l29l) obtained without the CS term. 
Rather, the appearance of knotted dual Abrikosov vortices is a qualitatively novel effect, 
which does not exist in the absence of a CS term. 

Notice flnally that the topological results obtained in this paper can be generalized to the 
cases of SU(3)- and SU(A^c)-inspired dual Landau-Ginzburg-type theories along the lines of 



Ref. 
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Appendix A: Some details of the derivation of Eq. (128p 

With the use of Eq. (l26l) . and owing to the conservation of j^, one has 
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+ if^vf [ [di [ D^,r>) (di [ orr, ) . (m) 

x,y \ J z / \ J u 



We furthermore assume the standard normahzation (1) = 1 of the functional average, which 
imphes a division by the functional integral j T>J^e~^ ■''^•y corresponding to the 

first term on the right-hand side of Eq. (lAip . Thus, we always imply that the measure PJ^ 
is normalized by a division by this integral. 

The last term in Eq. ( lAll) can be represented, through the integration by parts, as 
(/ir^)^ ^ D^j^ Dq^J^. The y-integration in this expression is straightforward: 

gip(x-y) r gjq(y-u) r gip(x-u) ^ 



y -J p 



^m'^o'" = / / ^—-2 / — = / 2f 2^ 2^ = " (^o" " , (A2) 

' 'p p2 + m2 q2 p2(p2 + 



where the equality p2(p2^_|_^2-) = ^ — ^r^^j has been used at the last step. Using 
further the explicit form of m^, Eq. ( l27l) . we can represent the last term in Eq. ( lAlll as 

In the second term on the right-hand side of Eq. ( lAll) . one can use the equality 
J^Dq^J^ = Dq^O^JI, which yields the same y-integration as in Eq. flA2p : J^D^Dq^ = 
(-Dq^ — D^). Upon the subsequent integration by parts, we obtain for this term the fol- 
lowing expression: ^e^^A y J^Ju^x {Dq^ — D^). Noticing also the definition of the Gauss' 
linking number, J) = Ef^^x J^^ J^HO^Dq^ , we arrive at Eq. ( l28l) . 

Appendix B: Some details of the derivation of Eq. (|39p 

For 0's obeying condition (157|) . one can use relations (ITS]) and ( 1201) to obtain the inequality 

« (Bl) 



02^ 

Owing to this inequality, the term ■^2^{d^j_ip^Y Eq. (138|) can be neglected in comparison 
with the absolute value of the term — g^^-^^^. The resulting Gaussian ^-integration can be 
performed by seeking the saddle-point function in the form ip^ = (f^}^ + i(p^^\ and solving 
the so-emerging system of equations for (fll'^ and f^^^. The result can be written as 



where = 2txJ^ + = -;^e^uxdvix are respectively the real and the imaginary parts 

of the current which couples to in Eq. (138|) . The obtained real and imaginary parts of 
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the saddle-point function Lp^ entering Eq. (lB2p read 



and 




with the new mass parameter M. = Furthermore, in the hmit (IBip at issue, the 

0{M. / m)-teims in Eqs. ( 1B3P and ( 1B4I) should be neglected compared to 1. That yields the 
following saddle-point expressions for y^^^'' and y^^^^ 

A'" J y J y 

Substituting them into Eq. (lB2p . we obtain for the Wilson loop in the limit (I37|) expres- 
sion ([2nD- 
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